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Abstract 

Based on two different iteration procedures the groundstate wave functions and 
energies for N-dimensional generalized Sombrero-shaped potentials are solved. Two 
kinds of trial functions for the iteration procedure are defined. The iterative so- 
lutions are convergent nicely to consistent results for different choices of iteration 
procedures and trial functions. 
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1. Introduction 



Recently the generalized radially symmetric Sombrero-shaped potential in 
A^-dimensional space is proposed by R. Jackiw[l]: 

V(r) = \g\r'-rl)\r' + Arl), (1) 

where = (2 + A^)/3, and A are arbitrary constants. He also chal- 
lenged to apply the iterative method developed by R Friedberg, T. D. Lee 
and W. Q. Zhao[2] to solving this problem. His question is properly an- 
swered by the three authors in Ref. [3] for one-dimensional case. The same 
problem will be solved for the A^-dimensional generalized Sombrero-shaped 
potential in this paper. The corresponding Schroedinger equation for the 
groundstate radial wave function is 

with k = {N — l)/2. The boundary conditions are 

V^(oo) = ^'(0) = 0. (3) 

When g = 1 and A = 2 the solution of the groundstate has an analytical 
form as ip{r) = with the eigenvalue Eq = Tq. However, for arbitrary 

g and A the groundstate wave function has no analytical form. In the 
following two iterative solutions for the groundstate of (2) are presented. To 
apply the iterative methods we introduce the trial function 0(r) satisfying 
another Schroedinger equation 

+ Vir) - h{r))ct>{r) = gE,cl>{r) = {E - A)0(r), (4) 

where h{r) and A are the corrections of the potential and the groundstate 
energy. Starting from this trial function 0(r) we perform two iterative 
procedures. Define the exact wave function as 

^(r) = /(r)0(r) = e'^^^'Uir). (5) 

The iteration performed for /(r) and A is named as /-iteration [2] and the 
one performed for r(r) and A is named as r-iteration[4] in this paper. For 
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the /-iteration two iterative series of {/n(^)} and {A^}, n = 0, 1, • • • are 
introduced with /o(r) = 1 and Aq = 0. For the r-iteration two iterative 
series of {T„(r)} and {A^}, n = 0, 1, • • • are defined with ro(r) = and 
Aq = 0. The iterations for these two sets can be performed according to 
the following equations. For /-iteration we have [2] 



OO 

/ r^V^(r)/i(r)/„_i(r)c?r 

= ' ^ . 7\ ' 





r J y 



fn{r) = Utc) - 2 / -J^^ J 2:''0'(^)(A{ - h{x))U-i{x)dx, (66) 

where rc could be chosen as rc = or rc = oo and the normalization is 
chosen as fn{fc) = 1- As for r-iteration we have[4] 



oo 



/ r^V'Cr) {h{r) - \ {T'^_^{r)f)dr 

= oo ^ , (6c) 

/ r^'^(l)^\r) dr 

t'^{t) = 2T-^H-\r) ]y'H\y)W^ - h{y)) + \{T'^_Mf\ dy M 

^ 

where T'(r) = ^. The detailed derivation is summarized in Appendix. To 
ensure the convergency of the iterative methods it is necessary to construct 
the trial function in such way that the perturbed potential h{r) is always 
positive (or negative) and finite everywhere. Specially, h{r) — ^ when 
r — ^ oo. In the following we construct two different trial functions for the 
iteration procedures. 



2. Trial Functions 

Trial Function I 

Introduce 

0(r) = e-^"^'"^ (7) 
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Now substituting (7) into (4) we obtain the equation for So{r): 

9k 

S'o(ry- - -S',{r) - S?(r) = 2(V(r) - h{r) - 9E0). (8) 

Therefore 

1 Q.k 

h(r) + gE, = V(r) - -(S',(rf - -S',(r) - SJ(r)). (9) 

For a finite h{r) it should not include terms with positive power of r. Since 
the highest order of r-power in the potential is 6 and V{r) has only even 
powers of r, we first assume 

Soir) = (ar"^ + cr^ + e) + mlog(ar^ + 1). (10) 

Substituting (10) into (8), to cancel term we have a = g/4. Coefficients 
e and a only change the normalization and we simply set e = 0, a = 1 and 
finally obtain 

S'o(r) = (|r^ + cr^)+mlog(r2 + l). (11) 
To cancel the terms with and we set 

c = ^giA - 2)rl and m = ^ig + 3)ri - ^giA + 2f4 (12) 



and obtain 
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h(r) = 2m(m-\-l)—^ — —7:-^(rng(A — 2)rQ — 2Tng-\-2rnk — 2Tn —m) 

(r^ + 1)^ + 1 

(13) 

gEo = \Aghl + 2mg + (2/c + 1 - 4.m)-^g{A - 2)rl. (14) 

When g = 1 and A = 2 we have m = 0, /i(r) = and the trial function 
is just the exact solution of the Schroedinger equation. To look in more 
details the behavior of the trial function we choose g = 1 and see the 
change of the trial function with the parameter A. When A = 2 our trial 
function is just the exact solution with the maximum of the wave function 
at r = 0. For A 2 we always have h{r) < and finite, and h{r) — ^ 
when r — ^ oo. This ensures the convergency of the iterative procedure. 



When A < 2 the potential is more centered at r = and the trial function 
keeps its maximum at r = 0. When A > 2 the potential is more like a 
double- well and the trial function has maxima at r 7^ 0. 



Trial Function II 

we can introduce another trial function 

(/>(r) = {j^fe-^^^^^^-^^^^^ (15) 

satisfying the Schroedinger equation (4) and the boundary condition 

0(00) = 0. 

The parameter a in (15) is fixed to ensure the boundary condition 

^'(0) = 0, 

namely 

^<sj(o) + <s;(o) + - = o. (16) 

a 

Substituting (15) into (4), we compare terms with the same power of g. 
From p^-terms we obtain 

S'oir) = = (r^ - rl)^r^ ^ Arl (17) 

To ensure the h{r) satisfying the convergence condition, Si{r) is defined 
in a special way to prevent terms with positive powers of r presenting in 
h{r). For terms we have 

--(2S'oS[-SS-2ka^/^^^TA4] - ka^ = Eq. (18) 

Z V / r=ro 

Introducing 

Eo = + + (19) 

and defining 

= -ka'^ (20) 

we write 

S'oS[ = i^S'^ - Ei'^) + (fca^^^^TA^I - eP). 



Since 5*0 (ro) = we obtain 

4'^ = Is'^iro) = rlVlTA, (21) 
E^^^ = karoVl + A (22) 

and 

S'l = il^'^ - 4')/^o + (fcaV^^M^ - Et)/S',. (23) 
Substituting «So(r) into (23) we have explicitly 

c/.x_ r^ + {l + A)rl 



^^2 + Ar§(rVr2 + Arg + r§^/^TA) 2(r2 + Arg) 

ka 



'r2 + ^r§(^r2 + Ar^ + roVl + ^) 
The expression for h{r) is 



(24) 



2^ 2 (r + a)2 r(r + a) ^ r(r + a) 



r(r + a) r(Vr2 + Ar^ + r) 

Substituting 5'o(0) and S'i{0) into (16) we obtain an equation for the pa- 
rameter a 

ka^ + iro VTTA -grlA){VA + VTTA) a + krliA + ^A{l + A)) = 0. (26) 

For the above equation to have real solutions of a the following restriction 
is put on the parameters g and A: 

{VTTA - gA4frl{^/A + VTTA) - A^Ak'^rl > 0. (27) 

For example, (27) requires g > 0.922 when A = 2 and A > 1.81 when 
g = 1. When (27) can not be fulfilled the condition 0'(O) = can be 
satisfied by introducing the trial function as 

4>revir) = (t){r) + ^(t)-{r) for r < rg (28a) 



and 

(Previr) = (1 + e0_(ro)/0(ro))0(r) for r > ro (286) 
where (j)-{r) is defined as 

0_(r) = {j^f e-^^^^-'^-^^^'\ (29) 

The parameter ^ is fixed to satisfy the condition 0'^e^(O) = 0, namely 

0'(O)+^0'_(O) = O. (30) 
Correspondingly the Schroedinger equation satisfied by 0rev(^) is 

+ ^(^) - ^re^(^))^re.(r) = y^O0re^(r), (31) 

where hrevif) = h{r) for r > ro and 
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V r[r + a) 

-ka'—j==^^—)4'-(r)l<t>rev(r) (32) 

for r < tq. It is interesting to notice that the conditions (26) and (30) for 
</>'(0) = also ensure h{r) and hrev{r) to be finite when r 0, which is 
necessary for the convergency of the iteration procedure. 
By integrating (17) and (25) we obtain *So(r) and Si{r) as 

SQ{r) = \r\lr'^ + ylr§(2r^ + Arl - 4ro) - \{A^rl + 4^ro) ln(r + ^r^ + ^rg) 
8 8 

(33) 

5,(r) =ln(r+r„)+-ln(r +Ar„)+(-+-) ln^==^_j— 

(34) 

Substituting them into (15) or (28)- (29) gives the final expression of the 
trial functions. From (25) we can also reach 

1/^/2 c^n^ J , ka i 



is — S"\ = - h 

2^ ^ 8{r'^ + A4y{a + (3) 2r{r + a){r'^ + Arl)a' + (3' 



fcV 

^2(r2 + Arl){^Jr^ + Arl + r^Vl^f 
kari2^r^ + Arg + rpVl + A) 

with 

7 = 225r^ + 270(1 + 2A)rV^ + 3(188A^ + 216A - 5)rVJ 

+36^(8^2 + 10^ - l)rV§ + 4^2(4^ + 
y = 9r^ + 3(4A + l)rV^ + 4A(1 + A)rJ, 
a = 15r^ + (18^ - 6)rV^ + (8^^ + 12A + 7)r2r(; + (SA^ + 2^)r§, 
= 8VlTArlr(3r^ + (2A - l)r^yr2 + Arg 

and 

a' = r(3r^ + (2A-l)rg), 

Substituting (35) into (25) and (32) gives the final expressions of h and 
hrev With above results for /i, and hrev, 4>rev we are ready to perform 
the iteration procedure. 

3. Numerical Result 

Starting from the above defined two sets of trial functions 0(r) and the 
related /i(r), we can perform the iterations based on /-iteration of (6a) 
and (6b) or r-iteration of (6c) and (6d). Our numerical results show that 
although the two iteration procedures look quite different and the two trial 
functions are defined in different ways the finally obtained wave functions 
and eigenvalues for the groundstate convergent nicely to the same final 
shapes and values. Now we give some more detailed discussions about our 
results. Let us take A/" = 3 as an example. 

For g = l and A = 2 
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For the trial function I, when g = 1 and A = 2, as mentioned before, 
the trial function gives just the exact solution of the groundstate 

with Eq = Tq. However, for the trial function II, it is necessary to fix the 
parameter a in the trial function (15) first by solving (26). It gives 

a = 4.4267 or a = 1.2976 

Performing the iteration based either on /-iteration or on r-iteration the fi- 
nal convergent result of the wave function and the eigenvalue of the ground- 
state is consistent to the exact solution. The trial function and the final 
exact wave function for the groundstate is plotted in Fig. 1. It is inter- 
esting to observe the transition of the shape of the wave function for the 
trial function with maxima at a finite r to the final convergent one with 
only one maximum at r = after the iteration procedure, as the exact 
groundstate wave function should be. This answered the question raised 
by R. Jackiw[l] in A'^-dimensional case: Even the trial function proposed 
has its maxima at r > the iteration procedure would still reach the exact 
solution of the groundstate wave function with its only maximum at r = 0. 

Comparison of the Two Iteration Procedures 

In Table 1 and Table 2 the eigenvalues of the groundstate obtained from 
the r- and /-iterations are listed respectively, based on the two different 
trial functions for different parameters g and A. Comparing the two iter- 
ation procedures, it can be seen that r-iteration is convergent faster than 
/-iteration. The numerical calculation takes also less time to reach the con- 
vergent result for r-iteration. This can be understood by comparing the 
formula (6a) and (6b) for /-iteration with (6c) and (6d) for r-iteration. 
First, in the formula for the energy correction, the denominator changes 
in each order in /-iteration while it needs only to calculate once for the 
whole r-iteration procedure. Besides, one fold less of integration is needed 
for each order of iteration in the r-iteration procedure since it is related 
only to r^. These two advantages speed up the numerical calculation of 
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r-iteration very much. Although the two iteration procedures look quite 
different with different convergent speed, it is shown clearly in the two ta- 
bles that the two iteration procedures do give the same convergent results. 

Comparison of the Two Trial Functions 

For the two trial functions, the trial function I is closer to the ground- 
state and needs less orders of iteration to reach the exact result in most 
cases. This can be seen clearly for the case of g = 1 and A = 2. The trial 
function I gives already the exact solution while for the trial function II 
which has maxima at r > the exact groundstate wave function with its 
only maximum at r = can be reached only after the iteration. In fact, for 
different parameters g and A the trial functions I always have shapes sim- 
ilar to the exact solution, while the trial functions II differ from the exact 
ones in their shapes for p < 1 and A < 2. Although the iteration process 
for the two trial functions is quite different, the two iteration procedures 
with the two sets of trial functions always reach the same final results of 
eigenvalues and groundstate wave functions. 

Change of the Wave Function Shapes with Parameters 

As examples the obtained groundstate wave functions after the iteration 
procedure are plotted in Figs. 2 and 3 for A = 2 and p = 0.5, 1 or 2, and 
for p = 1 and A = 1, 2 or 3, respectively. It is interesting to see the 
transition of the form of the obtained groundstate wave function from the 
shape with maximum at r = to the one with maxima at a finite r, 
becoming a degenerate groundstate, when g increases from < 1, passing 
1 to > 1 for A = 2, or when A increases from < 2, passing 2 to > 2 for 
g = 1- The results seem to show that the groundstate wave functions in 
the region g < 1 and A < 2 have the shape with only one maximum at 
r = 0, while in the region outside the wave functions become degenerate 
at a finite r. Their maxima move to larger r when the parameters g and 
A increase further. 
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Appendix 



For a particle with unit mass, moving in an N-dimensional potential V{r)^ 
the ground state wave function ^(r) satisfies the following Schroedinger 
equation: 

(_iv2 + V{r))^{r) = E^{r). {A.l) 

Introduce a potential correction h{v) and define another wave function $(r) 
satisfying the Schroedinger equation 

(_lv2 + F(r) - /i(r))$(r) = {E - A)$(r), (^.2) 

where A is the corresponding energy correction. Multiplying (A. 2) on the 
left by ^(r) and (A.l) by *l*(r), their difference gives 

_lv(W^^ - W$) = -{h - A)$^f. (A3) 

Let 

^^(r) = $(r)/(r) = ^{y)e-<''\ (A.4) 

In Ref.[2], an equation is introduced for (/(r), A). Here according to 
the same procedure similar results are deduced for (r(r), A). For the 
convenience of comparison we list both results together in the following. 
Substituting (A. 4) into (A. 3) equations for (/, A) and (r, A) could be 
obtained as following: 

-lv($M) = (-/^ + A)$V (A5) 

and 

iv($2Vr) = (-/i + A + ^(Vr)2)$2. (^.5) 

The integration over all space for the left hand side of above equations is 
zero. This leads to the expressions of the energy correction related to / 
and T, respectively: 
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and 

Jrfr$>-l(Vr)^] 

^ 7d7*5 ■ 

Introducing 

Al, A^2, a;, i = foT T (A9) 

and 

/i, /2, fn, ■■■ or n, r2, T„, ••• (^.10) 

the two iterative series are defined as 

iv($Vr.) = {-h + a; + ^(Vr„-l)2)$^ (^.11) 

- Jd^^ • ^^-^^^ 

and 

-^V($V/„) = {-h + A/)$Vn-i, (^.13) 

A-' = (A. 14) 

For later convenience the iteration series for {r„, A^} defined in equations 
(A. 11) and (A. 12), originally introduced in Ref.[4], is named as the r- 
iteration, while the one for {/„, A{} given in (A. 13) and (A. 14), originally 
introduced in Ref. [2], is named as /-iteration in this paper. By introducing 
the external electrostatic charge distributions 

< = + a; + l(Vr._i)^)$^ (A15) 

and 

a^ = (-/i + Ai)$Vn-i, (A16) 
correspondingly also defining displacement electric fields 

K = \^^Vt^ {A.17) 

and 

Di = -l^'^fn, {A.18) 
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where one can define k = $^ corresponding to tlie dielectric constant in 
the usual electrostatic problem. Now equations (A. 11) and (A. 13) can be 
expressed by the Maxwell equations of the electrostatic analog problem for 
r- and /-iteration [2, 4]: 



V-i^; = <, for i = r, /. 



(A19) 



If this electrostatic analog problem (A. 19) can be solved numerically, start- 
ing from the initial conditions 

A5 = 0, TO = 



or 



Ai( = 0, 



^0 — ^5 /o — 1, 

the corrections of the groundstate energy and wave function can be solved 
by iteration procedures (A. 19) and (A. 12) or (A. 14). However, the dielec- 
tric constant = in this problem is not a constant, but changes with r. 
This makes it quite complicated to solve (A. 19). Works along this direction 
are still in progress. 

For radially symmetric potential V{r) and potential correction h{r) the 
problem for solving the groundstate can be simplified and is related only 
to the radial variable r. By separating the angular variables[2] (A.l) and 
(A. 2) can be reduced to equations (2) and (4) for the groundstate radial 
wave functions 'ip{r) and 0(r). Multiplying (4) on the left by '0(r) and (2) 
by 0(r), their difference gives 



dr 



dr 



(^.20) 



Let 



ip{r) = (f){r)f{r) = 0(r)e 



{A.21) 



The equations for (/(r). A) and (r(r). A) are deduced as following 



d_ 
dr 



dr 



and 



d 
dr 



dr 



= 2r 



2k 



h-A- - — 

2\drJ 



(A.22) 



(A.23) 
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The integration of the left-hand side of Eqs. (A. 22) and (A. 23) over r = 
to oo is zero, which gives the expressions of the energy correction 

/ r^'^(f)^{r)h{r)f{r)dr 
A = , (A.24) 

/ r^^(f)^{r)f{r)dr 

and 

00 

/ r2V^(r) {h{r) - i (r'(r))2)dr 

A = 53 (^.25) 

/ r'^^(/)'^(r) dr 

with t' = ^. Introducing the two iterative series {/n(^)}, and 
{r„(r)}, {A^} with n = 0, 1, • • • it is easy to obtain the iteration equations 
(6a)-(6d) from (A.22)-(A.25). 

Acknowledgement 

The author would like to thank Professor T. D. Lee for his continuous 
guidance and instruction. 



References 

1. R. Jackiw, Private communication 

2. R. Friedberg, T. D. Lee and W. Q. Zhao, Ann. Phys. 321(2006)1981 
R. Friedberg and T. D. Lee, Ann. Phys. 316(2005)44 

3. R. Friedberg, T. D. Lee and W. Q. Zhao, arXiv: 0709.1997, 
Ann. Phys. (2007)doi:10.1016/j.aop.2007.09.006 

4. Zhao Wei-Qin, Commun. Theor. Phys. 43(2005)1009 



14 



Table 1. Eigenvalues of groundstates for A/" = 3 based on r-iteration 



A 

9 A 


irial mnc. 


T — 1 

^0 


7 — 1 

El 


T — 1 


T — 1 

^3 


7—1 

E4 


7 — 1 


0.5 2 


I 


1.1629 


1.3978 


1.3763 


1.3772 


1.3773 




0.5 2 


TT 

II 


-0.4300 


1.3963 


1.3763 


1.3773 


1.3773 




0.93 2 


I 


2.0237 


2.0352 


2.0351 


2.0351 


2.0351 




0.93 2 


TT 

II 


2.0921 


2.0457 


2.0355 


2.0351 


2.0351 




1 2 


I 


2.1517 












1 2 


II 


-(S.6479 


2.1523 


2.1517 


2.1517 


2.1517 




2 2 


I 


3.6066 


4.1140 


4.1093 


4.1094 


4.1094 




2 2 


TT 

II 


5.5581 


4.1362 


4.1123 


4.1097 


4.1094 




1 1 


I 


2.5073 


1.8400 


1.8392 


1.8392 


1.8392 




1 1 


II 


-2.3537 


1.8920 


1.8330 


1.8394 


1.8393 


1.8392 


1 1.9 


I 


2.1225 


2.1215 


2.1215 


2.1215 


2.1215 




1 1.9 


II 


-3.8095 


2.1232 


2.1215 


2.1215 


2.1215 




1 3 


I 


3.5310 


2.4630 


2.4426 


2.4418 


2.4418 




1 3 


II 


3.6773 


2.4675 


2.4437 


2.4419 


2.4418 





15 



Table 2. Eigenvalues of groundstates for A/" = 3 based on /-iteration 



A 

9 A 


irial mnc. 


T — 1 

^0 


7 — 1 

El 


T — 1 


T — 1 


7 — 1 


7 — 1 


0.5 2 


I 


1.1629 


1.3978 


1.3705 


1.3786 


1.3770 


1.3773 


0.5 2 


TT 

II 


-0.4300 


1.3963 


1.3795 


1.3775 


1.3773 


1.3773 


0.93 2 


I 


2.0237 


2.0352 


2.0351 


2.0351 


2.0351 




0.93 2 


TT 

II 


2.0921 


2.0457 


2.0337 


2.0352 


2.0351 


2.0351 


1 2 


I 


2.1517 












1 2 


II 


-(S.6479 


2.1523 


2.1516 


2.1517 


2.1517 


2.1517 


2 2 


I 


3.6066 


4.1140 


4.1088 


4.1094 


4.1094 




2 2 


TT 

II 


5.5581 


4.1362 


4.0976 


A -1-1 /^O 

4.1108 


A -1 f\r\c\ 

4.1092 


A -1 A 

4.1094 


1 1 


I 


2.5073 


1.8400 


1.8392 


1.8392 


1.8392 




1 1 


II 


-2.3537 


1.8920 


1.8473 


1.8402 


1.8393 


1.8392 


1 1.9 


I 


2.1225 


2.1215 


2.1215 


2.1215 


2.1215 




1 1.9 


II 


-3.8095 


2.1232 


2.1214 


2.1215 


2.1215 


2.1215 


1 3 


I 


3.5310 


2.4630 


2.4464 


2.4425 


2.4419 


2.4418 


1 3 


II 


3.6773 


2.4675 


2.4353 


2.4425 


2.4417 


2.4418 
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Figure Caption 



Fig. 1 Trial Function 0(r) and Groundstate Wave Function ip{r) 
for AT = 3, p = 1 and A = 2. 

Fig. 2 Groundstate Wave Function '0(r) for N = 3, g = 1 
and A = 1 (thin) , 2 (middle) and 3 (thick) . 

Fig. 3 Groundstate Wave Function ip{r) for N = 3, A = 2 
and <7 = 0.5 (thin), 1 (middle) and 2 (thick). 




Fig. 1 Trial Function 0(r) and Groundstate Wave Function ^(r) 

for = 3, ^ = 1 and A = 2. 



17 




Fig. 2 Groundstate Wave Function ip{r) for N = 3, g 
and A = 1 (thin), 2 (middle) and 3 (thick). 




Fig. 3 Groundstate Wave Function ^(r) for N = 3, A 
and g = 0.5 (thin), 1 (middle) and 2 (thick). 
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